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Abstract. Conrey, Farmer, Keating, Rubinstein and Snaith have given a 
recipe that conjccturally produces, among others, the full moment polynomial 
for the Riemann zeta function. The leading term of this polynomial is given 
as a product of a factor explained by arithmetic and a factor explained by 
combinatorics (or. alternatively, random matrices). We explain how the lower 
order terms arise, and clarify the dependency of each factor on the exponent 
k that is considered. 

We use extensively the theory of symmetric functions and representations of 
symmetric groups, ideas of Lascoux on manipulations of alphabets, and a key 
lemma, due in a basic version to Bump and Gamburd. Our main result ends 
up involving dimensions of skew partitions, as studied by Olshanski, Regev, 
Vcrshik, Ivanov and others. 

In this article, we also lay the groundwork for later unification of the combi- 
natorial computations for lower order terms in the moments conjectures across 
families of L-functions of unitary, orthogonal and symplcctic types. 



1. Introduction 



1.1. History. This paper concerns the Riemann zeta function and its moments. 

In 2000, Keating and Snaith [KSOOa, KSOOb computed moments of characteristic 
polynomials of matrices in the unitary group U (N) and suggested the following 
conjecture. 

Conjecture 1 ((see |KS00b| )). For any positive integer k, 
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^This is usually referred as the random matrix theory factor. We feel however that the adjective 
"combinatorial" is more appropriate, particularly in light of the results presented here. 
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where A g is the characteristic polynomial of a matrix g distributed under Haar 
measure. 

The normalization by (fc 2 )! is there to guarantee that gk £ N. 

Their conjecture followed a few unconditional results early in the 20 th century 
(Hardy and Littlewood for k = 1, Ingham for k = 2) and a few conjectures in the 
late 20 th century (folklore for the exponent of log T and the general shape of the 
coefficient, Conrey and Ghosh |CG84j for more conjectural information about this 
coefficient, Conrey and Ghosh and Conrey and Gonek CG98, CG01 for the conjec- 
tured values of (73 and 34). While agreeing with all the previously known results or 
conjectures, Keating and Snaith's contribution was to identify the sequence gk that 
should work: (&) = (1, 1, 2, 42, 24024, 701149020, 1671643033734960, ■ ■ ■ ). This is 
sequence A039622 in [A0311J . 

Soundararajan has shown that, assuming the Ricmann Hypothesis, the exponent 
of logT in Equation (jTJ) is correct, at least up to an e (see |Sou09j for the precise 
statement). 

The work of Keating and Snaith attracted significant interest, and led to several 
simplifications or new interpretations. 

In 2002, the collaboration of Conrey, Farmer, Keating, Rubinstein and Snaith 
[CFK+05llCFK+08] gave a recipe leading to a sharpening of Conjecture ((l) for inte- 
gral k. This is based on earlier unconditional work of Motohashi |IM06j sharpening 
the theorem of Ingham for k — 2. Their final conjecture takes the form 

Conjecture 2 ((see [CFK+05] . (CFK+08) )). Let k e N. There exists a sequence 
of polynomials, of degree k 2 and with leading coefficient equal to jpfy , given either 
im plicitly by a Ik-fold contour integral in [CFK + 05] or explicitly by Theorem 1.2 
in [CFK+08] . such that, for all e > 0, 

2k 

(A) 







Pk ( log ^ ) 
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The information on the leading term shows that Conjecture [2] implies Conjec- 
ture [TJ While the leading term is simply given as a product of the ak and the gk 
factor, it became a major issue to understand the interaction of the arithmetic and 
combinatorial terms for the lower order coefficient^. 

Diaconu, Goldfeld and Hoffstcin |DGH03j immediately reformulated these con- 
jectures (or their leading order) as consequences of further conjectures in the then- 
emerging theory of Multiple Dirichlet Series. This is in some ways even closer to 
the works of Motohashi. 

In 2003, Beineke and Bump |BB04a[ IBB04bj remarked that the CFKRS recipe 
could potentially be explained as constant terms of Einsenstein series on GL(2fc), 
where a distinguished role for the subgroup GL(fc) x GL(fc) is to be expected. 

In 2004, Bump and Gamburd [BG06j simplified Keating and Snaith's computa- 
tions on moments of random matrices. This work stressed the combinatorial nature 
of gk, which has no arithmetical content. 

1.2. Motivation. Conjecture [2] is the most precise conjecture on moments of the 
Riemann zeta function that has been formulated so far (and is part of a large set of 
other similar conjectures on families of L-functions) . The polynomials Pk are thus 



2 As far as is currently known, there is no expression for these lower order coefficients directly 
in terms of random matrices, thereby justifying the footnote on page 111 
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very interesting, but they need to be known very precisely to be useful. Indeed, they 
tend to have small leading coefficient and comparatively large middle coefficients, 
which is problematic for testing: for actual numerical data (i.e. large but finite T 
[HOIOj ). it is harder to tell which term would contribute most to P/ £ (logT). 

The heuristics leading to Conjecture [2] are based on approximations of zeta by 
Dirichlct polynomials, followed by complicated combinatorial manipulations. This 
is similar to the earlier conjectures of Conrey and his collaborators. Unfortunately, 
the end result in [CFK+05] for Pk is extremely implicit, given as the residue of a 2k- 
fold integral, which can be either obtained symbolically or integrated numerically. 
In either case, limitations creep in that drastically bound the size of achievable k. 
This formula (a 2fc-fold integral) makes it also very difficult to think of the analytic 
continuation in k of the Pk , which should exist (these questions are tied to bounds 
for the Riemann zeta function on the critical line, and eventually to the Lindelof 
Hypothesis). The end result in [CFK+08] answers some of these questions, but 
the computations are even more complicated and in the end useful numerically but 
not so illuminating. In particular, a polynomiality result for components of the 
Pk is proved in their Theorem 1.3. These polynomial components, which they call 
iVfc(a;/3), are then recovered using either interpolation or dctcrminantal formulas. 
Concern for these polynomial components end up occupying a significant part of 
[CFK+08) . but in fact we would argue that the Nk(a; /3) are not satisfactorily 
explained by their paper. As evidence we will point to a very recenlQ preprint 



|GHRR12] which aims to extend the treatment of |CFK + 08 to symplectic and 
orthogonal families. While the arithmetic side requires little new ideas, the study 
of what they now call N\(k) (in their Theorem 1.1) and again its polynomiality is 
the hardest part and requires to perform complicated calculations anew. 

Starting with the recipe of jCFK + 05] . we wish to provide an alternative to the 
computations of [CFK+08] . avoiding entirely the use of multiple integrals. 

1.3. Statement of the main result. In accordance with CFKRS, we define func- 
tions c r (k) at integers k,r > such that 

(5) Pk(x) =c (k)x k2 +c 1 (k)x k2 ~ 1 H hcy(fc). 

Before we give a formal statement of our main results, we qualitatively describe 
them. The narrow view of our results is that the polynomial components (or in fact 
linear combinations of the Nk(a; (3)) have a combinatorial interpretation counting 
standard tableaux of skew shape or alternatively paths in the Young lattice. This 
interpretation, which is very natural to a combinatorialist, immediately leads to 
the polynomiality results thanks to the theory of Frobenius-Schur functions, and 
gives numerous efficient dctcrminantal formulas to compute them. This interpre- 
tation will even give the analytic continuations in k. The broader view is that this 
interpretation will turn out to work in the symplectic and orthogonal cases (in the 
upcoming paper Combinatorics of the lower order terms in the moment conjectures 
for the symplectic and orthogonal families of L -functions), thereby unifying our 
understanding of the computations needed for lower order terms across the three 
types of families. 

We comment further on the significance of these results after stating them for- 
mally. 



^In fact, a first version of this paper (with minor differences except in the introduction and 
Section 18.21 1 preceded their preprint on the arXiv. 
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Theorem 3. The coefficient CAr(fc) satisfies the equation 

(6) Cjv(fc) = ^ fc2 - jjy d K Adim(A,5 fe (K)), 

|«| + |A|=JV 

with the sum taken over pairs of partitions k, A of combined weight N , and 

• dim(/z, v), a nonnegative integer, the number of standard tableaux of skew 
shape \i \ v ( defined in Section \2.1.2\) : 

• Sk(n) a partition of size k 2 — \k\ (defined in Section^; 

• d K \ = d K \(k) a coefficient of arithmetic significance. 

To be more explicit, the d K \ are given in Proposition [?| Section \5\ as simple 
algebraic expressions (with rational coefficients depending on character tables of 
symmetric groups of size N, see the proof of Proposition^ in the W^, which are 
defined in Equation (|48p as 

r>l 

with the easily computable polynomials in k ( see Equation (|53[) ), of degree 
bounded by 2r, and the n th Taylor coefficient of the prime zeta function at r 
(see Section \2.4-2\ ). Moreover, the algebraic expression for the d K \ only involves 
Wjj,„ of smaller total weight in p, and v than N . 

The (summary of the) proof of this theorem is on page [1^1 where references are 
given to the different lemmas and propositions needed. Some of the main formulas 
are listed in the appendix, page 1251 

As will be clear from the definition of dim(^, v) and S^k), the formula would 
give cjy(k) = if applied when N > k 2 . 

The goal of this theorem is to express the cjv(fc) in terms of several functions that 
are very simple to precompute, with a clear sense of breaking down each into simpler 
and simpler functions and with many of them independent of k. We know from 
previous work of CFKRS that some components of the Pk will not be computable 
exactly, and we have relegated that part to computing the Taylor coefficients of the 
prime zeta function, which can be done to high numerical precision. The sum giving 

(n) 

the Wfjtu is the last bit that is problematic, as it is infinite. Note however that F r 
decreases exponentially fast in n to 0, so the tail of the sum could be bounded, for 
instance if we only want a numerical approximation. See also Section [9] for more 
details. 

Before commenting further on this theorem, we state another theorem that will 
clarify the dependency in k of the RHS of Equation (|6|) . 

Theorem 4. // k and A are partitions, then 

(7) dim(X,S k (K)) =dim(« > 5„(A)) = fc2 ■ (fc2 _ 1} ■ ggA ^ - ^ + ^ 

where B(k) is an explicitly computable polynomial in k of degree bounded by 2(|k| + 
|A|) and g], is the combinatorial factor in Conjecture^ 

The proof of this theorem is on page |2T)1 

We now return to the discussion of Theorems [3] and 0] in the context of existing 
literature. Theorem [3] has a structure similar to Theorem 1.2 of jCFK + 08| : both 
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give an explicit formula for the coefficient cjv(fc) as a sum over pairs of partitions 
of total weight N. In both theorems, each summand is a product of a term "ex- 
plained by arithmetic" and another one "explained by combinatorics" that has a 
nice dependency in k. The polynomials Nk(a,/3) in [CFK + 08] entered on the com- 
binatorial side, and this paper uses instead the function dim(A, Sk(n)). This has 
several advantages: 

(1) The paper jCFK + 08] offered to compute the polynomials Nk(a,j3) in k 
using either interpolation on the values of k x k determinants or via a 
more direct evaluation (in Equation (3. 1 1^0 of its Theorem 3.2) as a s x 
s determinant, where s can be bounded, after careful reading, by s < 

= with N the index of the coefficient CAr(fc) that is considered 
(remember from Equation (|5|) that it is meaningful to let N go up to fc 2 , 
when k G N). 

Since we rely on existing theory of the Frobenius-Schur functions, we 
can hope for better formulas to compute the dim(A, Sfc(«)). Indeed, the 
Giambclli formula in Equation (|55[) will give the dim(A, ^(k)) as determi- 
nants of size at most [V^VJ x |_"\/iV J , a significant gain from a computational 
perspective (see also Section 15721 for more on this bound). 

(2) Our formula gives a concise and direct formula for the analytic continuation 
in k (this requires Theorem [4] as well). 

(3) The derivation on the combinatorial side will be unified in Combinatorics 
of the lower order terms in the moment conjectures for the symplectic and 
orthogonal families of L-functions with the derivation for other families of 
L- functions. 

(4) On the arithmetic side, while [CFK+08] gives a succession of steps, we have 
figured out the exact final formulas. Numerically, this is an advantage: the 
non-exact part of the computation is limited to the evaluation of the Taylor 
coefficients iv , and then the truncation of the sum for W^ v . It should be 
easier to achieve better precision, and would be conceivable to use interval 
arithmetic to bound exactly the error actually made. 

(5) In [RHllj . Hiary and Rubinstein obtained uniform bounds on the coef- 
ficients cj\r(fc). For this, they had to bound the arithmetic side and the 
combinatorial side separately. We did no work in this direction, but com- 
ment on connections in Section 18.21 

1.4. Structure of the paper. This paper requires a variety of very different tools. 
In Section[2j we present the background needed on partitions, symmetric functions, 
alphabet manipulations and power series expansions of the Riemann zeta function 
and the prime zeta function. In Section [31 we present the beginning of the CFKRS 
recipe in details, but formulated using symmetric functions. In Section 21 we intro- 
duce the substitute to the 2fc-fold integrals of CFKRS, a lemma due to Bump and 
Gamburd in its basic form. In Section [5j we compute the final expression handed 
to us by the recipe, using that lemma. Section [6] contains the manipulations that 
allow us to transfer Dirichlet polynomials computations in the language of symmet- 
ric functions. In Section [3 we present the proof of Theorem HI Section |8j contains 
the various corollaries to the main results, and in particular the rederivation of 



Unfortunately, this formula has a typo on the right hand side: one needs to substitute s X s 
for k X k. 
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the Keating-Snaith leading coefficient. In Section [9j we explain how to compute 
the constants involved in our derivation. Due to the complexity of the formulas 
obtained, we have decided to add a table of the main formulas and definitions in 
the appendix, page 1251 

The author wishes to thank many people for useful discussions during the (long) 
period that led to this paper and in particular Profs. Bump, Farmer and Nikegh- 
bali. In addition, Profs. Conrey and Kowalski read preliminary drafts and provided 
feedback and corrections. Many tests and computations were performed using the 
computer algebra package sage jS + ll) . 

2. Some Essentials 

2.1. Essentials in the theory of partitions. We rely almost exclusively on def- 
initions and objects introduced in the first chapter of [M ac95j . We briefly present 
these definitions. 

2.1.1. Definition and easy invariants. An (integer) partition A = (Ai, • • • , A„) is a 
finite weakly decreasing sequence of non-negative integers, called parts. We define 
the weight |A| of A to be the sum ■ Xj. If this weight is n, we also use the notation 
Ahn. The length Z(A) of A is the maximal I such that A; ^ 0. An empty sequence 
() is thus also a partition, of length and weight 0. We denote this empty partition 

by 4>- 

We freely think of partitions in terms of Ferrers (also called Young) diagrams. 
See Figure [TJ page [15] for some examples. The Ferrers diagram of a partition is 
constituted of boxes. For each box □, we call the hook of □ the set consisting of 
the boxes exactly to the right or exactly below □, and □ itself. The hook length of 
□ is the cardinality of that set. 

The conjugate of A is the partition obtained by flipping the diagram of A along 
the main diagonal. It is denoted A*. 

Given a partition A, we can naturally associate a conjugacy class of permutations 
in S\\\ by considering that the Ai give the lengths of the cycles of permutations. We 
define z\ to be the order of the centralizer of an element of that conjugacy class. 
It is immediate to compute 



where Mj(A) counts the number of parts of length j (hence almost all of those will 



We also need what we call vectorized partitions. Given a partition A of length 
I < L, we sometimes want to consider the partition as a vector of length L. In 
that case we use the symbol A for the vector (Ai,--- , A;,0, ••• ,0), i.e. we pad 
the sequence of A with L — I zeroes. We say that A is vectorized to length L. 
We can of course add any two partitions vectorized to the same length: we do 
this componentwise, as vectors. We let pL be the partition (L — 1, • • ■ ,1) (with 
triangular Ferrers diagram), vectorized to length L (by appending one 0). 

In addition to conjugacy classes in symmetric groups, partitions are used to 
index many objects of combinatorial interest, such as symmetric function bases or 
characters of symmetric groups. We will explain some of this in Section l2~2l 



(8) 




be 0). 
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2.1.2. The dimension of a (skew) partition. We define the Young lattice as a graph 
(V,E). Let V n be the set of Ferrers diagrams of weight n. Set V to be Ll n >oV n . 
The (undirected) edges E only connect Ferrers diagram that differ by the addition 
or removal of exactly one box. Edges can thus only link partitions whose weight 
differ by 1. 

Given two partitions k, A, of weight k and I respectively, we let dim(K, A) to be 
the number of paths in the Young lattice of length exactly I — k connecting k and A. 
In particular, this will be if k > I, or if the Ferrers diagram of n is not contained 
inside the Ferrers diagram of A. 

There exist convenient formulas to compute dim(K,A), and these will be pre- 
sented in Section EH 



2.2. Essentials in symmetric function theory. We present here the concepts 
of symmetric functions. We follow entirely the conventions of Macdonald |Mac95| . 
except that we also define the symmetric polynomials po- 

2.2.1. The ring of symmetric polynomials. Let X = {xi,--- ,x n }. Consider the 
ring C[[xi, • • • , x n ]] of power series in the variables X. The symmetric group S n 
acts on this ring by permuting variables, and a power series is symmetric if it is 
invariant under this action. The symmetric power series form a subring 

(9) A n = A n (X)=C[[x 1 ,--- ,x„]] 5 ". 

We also set A„ = A„ n C[Af]. The ring A„ is a graded ring: we have 

(10) A„= ©A«, 

fe>0 

with each component consisting of polynomials of degree k (together with the 
polynomial) . 

We will use most often elements of A„, i.e. not just polynomials but infinite 
series of elements of A n . We call these symmetric functions. 

2.2.2. Bases. Combinatorialists have studied several different bases of polynomials 
in A n . We introduce them now. 

For each a = (ai, • • • , a n ) € N™, we set x a := x" 1 ■ ■ ■ x" n . Given a partition \x 
such that Z(/i) < n, we let the monomial symmetric polynomials be 

(11) nv(*) := £ x a , 

where the sum runs over all the distinct permutations a of A = (Ai, ■ • • , A n ). The 
set of xn\ such that l{\) < n and |A| = k forms a Z— basis of An\ 
For each r > 0, we define the complete symmetric polynomial f) r as 



(12) f) r :=^m A . 



Ahr 



In particular, we have l)o = 1. We also set f) M := Il'=i bi particular, we have 

For each r > 0, we also define the power symmetric polynomial p r as 

(13) Pr (X) := £ x r . 

xex 
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In particular^, we have po(^) = We also set p M := Jli=i I 11 particular, 
w^] have p,/, = I. 

The last set of polynomials is the Schur polynomials s M for fj, running through 
partitions. If l(fj,) > n, we set s fl (X) = and otherwise 

(14) 



after vectorizing fj, to length n if needed. Note that Sa = 1. 

2.2.3. Transitions between bases. For each pair of basis of A, there exists a transi- 
tion matrix from one to the other. They have long been known, and Macdonald 
devotes a whole chapter to them |Mac951 Chapter 6] . 

For each partition A there exists a character \ X °f the symmetric group S\\\. 
These characters are class functions, and thus only dependent on the cycle type of 
the conjugacy class considered. Denote the value of that character on the conjugacy 
class of cycle type given by the partition /i as x A (m)- This coefficient is known to 
be an integer, and the dimension of the character x\ *- e - X A (id) = X A ((1, ■ ■ ■ , 1)), 
equals dim(0, A). The whole theory is beautifully introduced in |Sag91| . 

It turns out these character values also serve as transition coefficients between 
the and the Sa- Indeed, we haveQ 

(15) s A = Y xX{ji)Ptl 
and 

(16) P P = X X Msx- 

AH Ml 

The only other transition we will need (actually a special case of the first transition 
presented above) is 

(17) bn(A)=J2-p K (A)- 

i Zk 

2.2.4. Pieri rule. This rule gives a multiplication formula in A for a Schur polyno- 
mial s K and the power polynomial pi: 

s K pi = X! Sk+ ' 

where n + runs through partitions \i such that dim(«, \x) equals 1. By repeated 
application of this rule, we see that 

s K (Pi) r = X! dim ( K >/4v 

/ih|«|+r 



^This is true even if some Sj is 0! Actually, we always avoid directly setting any Xi to 0, but 
we let them tend to 0. The value of p r {X) is then less shocking. 

^Macdonald does not define po, possibly because of the potential confusion arising from po ^ 
p^,. These will actually turn out to be quite helpful to us: they allow to temporarily hide away 
the dependency of expressions on the size of the set X . 

^Remember (from Equation (£Sj> ) that z M is the order of the centralizer of any permutation of 
cycle type fi, and that therefore the conjugacy class of cycle type fj, is of order 
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2.3. Essentials on alphabets. Wc only use a very elementary version of Lascoux' 
theory of alphabets. This paper could almost surely be simplified by using more of 
that theory. 

2.3.1. Definition. To us, alphabets will be finite sets of indeterminates. Sometimes, 
indeterminatcs might happen to take the same value, but they will always be dif- 
ferent indeterminates. If T> is an alphabet, its cardinality is thus constant, and 
denoted by \V\. 

We distinguish some special alphabets: 1 := {1} and := {0}, both of cardinality 

1. 

2.3.2. Notational shortcut for alphabets. Wc will use alphabets constructed on three 
main alphabets: A, B and T>. For simplicity of notation, we index elements of A 
by a, for B by j3 and for T> by 8. In the same spirit, we write 

n/(*) f ° r 

<5 SeT> 

Note also that there will be significant advantage in using the notation J\ s f(5) 
instead of = JXEi f^i) ^ or ms tancc: the latter notation imposes an (unnatural) 
order on the elements of T>, while the former localizes any dependency in the size 
of the alphabet, and spares the introduction of dummy indices. 

2.3.3. Basic operations on alphabets. Given an alphabet A, or maybe a pair of 
alphabets (A, B), one can define another alphabet X by performing some operation 
on the elements of the starting alphabet(s). We present a few constructions. 
Repetition. Given an alphabet A and k G N, wc define k ■ A by taking the union of 
k disjoint copies of A. By extension, for any k € C, we can define a new "abstract" 
alphabet by 

(18) p m (k-A):=kp m (A). 

Pairwisc addition. Given two alphabets A — {a} and B = {/3}, we define A®B := 
{a + 13 : a <G A, /3 E B}. This implies 

(19) p m (A®B) = J2 ()PMp m - u (B). 

Wc have for instance p Q (A®B) = p {A)p (B) = \A\ \B\ and p x (A © B) = \A\pi{B)+ 
\B\pi(A). Note how po turns out to be useful to obtain uniform formulas. 

This definition already includes the notion of translates of alphabets. Indeed, 
given an alphabet A — {a} and an element r £ C, we have A © {r} = {a + r}. 
Scaling. We can also scale alphabets by a complex number. We write this <fi r (A) := 
{ra : a € .4}, which implies 

(20) PmiM-A)) =r m p m (A). 

A special case occurs when r = — 1, where we use the shortcut A := <fi-i(A). 
We can also define scalings by a partition n: 



(21) 



<f> K (A) := (f> Kl (A) © • • • © Ki(K) (A). 
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Despite these simple definitions, it already becomes harder (but will be useful) 
to compute p m (^(i)). Let I = 1(k). Then, 

p m {<j> K (A)) = p m (<M-4) e • • • e <t> Kl(K) (A)) 



»l,-"iH 

m 



'i i ■ ■ j '(/ 



/ihm 



(22) = E (?) m M(«)M-A)W 



where f m ) := f m ) is a multinomial coefficient. The first equality follows 
from the definition of </> K (*4), the second from iterating the definition of pairwise 
addition, the third from the definition of pi(</) r (A)) and the fourth from the defini- 
tion of tn^(re) and re-sorting of the vector i into the partition /it. The fifth uses the 
definition of po(A). 

It follows immediately from this that we have Po(4>k{A)) = \A\ 1 ^ . 
Pairwise multiplication. Similarly to pairwise addition, we define pairwise multipli- 
cation by 

(23) A®B:={a-0:aeA,0€B}, 
which implies 

(24) p m {A®B) =p m (A)p m (B). 

Of course, scaling is a special case of pairwise multiplication. Indeed, 4> r {A) = 
{r} <E> A. In general, we have p Q (A ®B) = p (A)p (B) = \A\ ■ \B\. 
Elementwise exponentiation. Let l£R,X>0. Given an alphabet A = {a}, we 
set X A := {X a }. We immediately get: 

(25) Pm (X A ) := E ^^ PuiA). 

u>0 

This verifies p (X A ) = p a (A) = \A\. 
The definitions guarantee that 

(26) p m (X A ® B ) = p m (X A ® X B ) = p m (X A ) ■ p m (X B ) . 
Resolvent. Finally, we define 

(27) A(A;B):=l[(a-f3)= [J 7 , 
which is sometimes called the resolvent. 



2.4. Essentials in number theory. 
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2.4.1. The Riemann zeta function. For completeness, we redefine the Riemann zeta 
function 



n>l 



We arc actually interested in its power series expansion around 1, once the polar 
part is removed: 



C(i 



»)=-+ E^y 1 ^ 71 = - 1 + £(-ir v +1 

n=0 \ n=0 / 



where 7,; are the Stieltjes constants and thus 70 is the Euler constant. 
We will also need the expansion 

(28) \og( S -((l + S )) = -Y, [ —r9 C nS n , 

71=1 

with the Stieltjes cumulants, as defined by Voros in jVorlOl p. 25]. These co- 
efficients are more convenient for us, and are also more directly linked with the 
Riemann zeroes than the coefficients 7„ f |VorlQ[ p. 70] and }Vor031 p. 676]). 

The zeta function satisfies a functional equation £(s) = x(s)C(l~ s), with x( s ) = 
2 s 7r s ^ 1 sin (^) r(l — s) and more importantly to us 

(29) x{s )=(±Y + 



2-k J \ \t 

for s in a vertical strip in the complex plane. This leads to the relation 

-k / 1 \ / t \ -pi(B) 

2W 



(so) xQ+u) nx{\+«+p) = (lY 

|B|=fc 

also for s in a vertical strip. 



2.4.2. The prime zeta function. The prime zeta function is much less studied than 
the Riemann zeta function. Its definition is 

i) = V^ k s >i, 

which admits a meromorphic continuation to any simply connected domain sand- 
wiched between {s : SRs > 1} C Q C {.s : 5Rs > 0} and such that 51 does not contain 
any x with the pole or a zero of for a squarcfrec n |Fro68j . In any case, we 

o 

will only need C(s) for 5is > 1, at worst approaching 1. 

The prime zeta function is much more adapted to the computations that will 
follow. In fact, we need its Taylor coefficients at various integers. When r > 1, let 

(31) C('- + ^)=:E^ ( " )s ™' 

71>0 



withF r w =C(r). 
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When r = 1, we would like to do the same, but £(s) is not meromorphic at 1. 
The singularity can easily be characterized. This uses a Mobius inversion |Fro68j : 

/x(n) 
n 



(32) l(s) = J2^ log CM »*>1, 



n>l 

which seems to have already been known to Glaisher (1891). In a neighborhood of 
1, we can use this to study the singularity: 

(33) C(l + s) = log( S ) - £ ^Tf-^s" + E — l0 S C(n + ns). 

* — ' n! z — ' n 

n=l n>2 

This leads us to define the coefficients 

(34) C(l + s)=:-log(s) + E i? i" )s "' 

n>0 

with the observation that F^ = £„ >2 M™) 122 !^ = B i ~7o(~ -0.315718452), 
with J3i the Mertens constant. 

3. Following the CFKRS Recipe, at least initially 
CFRKS start with (an expression similar to) 

(35) ri c G +tt )r d * = ^). 5 f x G +tf ) VG + H d *' 

" — v — ' 

(*) 

where the functional equation is £(s) = x( s )C(l — s ) an d 2? — > (2k) ■ means that 
we are letting each of the 2k indeterminates of T> go to 0. 

The recipe is composed of three main steps. As CFKRS insist {CFK+051 p. 53], 
these are not individually heuristics: none of those steps can be properly justified, 
and some in fact give diverging expressions. Some large terms are added while some 
others are then thrown out. However, the conjecture is that the composition of the 
three steps gives the correct final formula. Indeed, it agrees with known theorems 
for special cases, and is otherwise reasonable and matches experimental data. We 
name the three steps as follows: 

(1) approximate functional equation step; 

(2) y-selection step; 

(3) diagonal selection step. 

We will always indicate a traditional equality with the traditional = sign. However, 
we also use the symbol =>■ to indicate one of the three unjustifiable steps above. 

We first state the approximate functional equation |Tit86| . Let < 3?s = cr < 1 , 
and t > C > with 2nMN = t. Then, uniformly in cr, 

N M 

(36) CM = H-s+ X« E ^ + O(error). 

n—1 n—1 

The approximate functional equation step replaces each C (h + 8 + it) factor 
in the RHS of (J35J| by 

1 /l , ..\r^ 1 



(37) 



E 
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which is cause for serious concern: neither of the two sums converge when the shift 
S is too small. Nevertheless, CFKRS carry on and replace each factor in (*) by 

433). 

Explicitly, wc obtain for (*): 



(*) 



(38) 



it 



n<a 



-k 

n 



6 + it dt 



E 



i 



7) 2 



+<5+it 



+ x 



dt. 



A naive expansion of this last product over S leads to a sum of 2 2fc terms. The 
alphabets A and B serve to encode which of the two terms is chosen for each factor. 
The sum is thus over the set partitions of T> = {5}. 



Eqn. ([38 



E IIE i+a+itTI 



A,B s.t. a n>l 
AuB=V 

AnB={} 



n>l 



p-it 



We now use information about \ contained in Equation (|29[) . For high value of t, 
on the critical strip, x(t) will oscillate very rapidly (due to the t lt term in Equa- 
tion (|2"9l ). This will lead to a negligible contribution once integrated over t. There- 
fore, the ^-selection step imposes the further restriction that \B\ = k. Indeed, 
those terms will dominate, due to the simplification of Equation (|30|) : 



(39) 




-Pi(B) 



A.B s.t. 
AUB=V 
|-4| = |B|=fc 



n e i+ Q +it ■ n e i- p-u f 

a n>l n p TO >1 J 



dt, 



where we have used the dummy variables m and n to prepare for what follows. We 
now expand the products over a and in the most naive way possible. This leads 
to the choice of m, ■ ■ ■ , associated to a±, ■ ■ ■ , ajt and mi, • • • , associated to 
\3\ , • ■ ■ , Pk- Each summand in the full expansion will be of the form 

i/ 



(40) 



c(n, m, a, /3; t) 



n 



m 



P, 



1 



n: 




It is immediate that 



(41) 



c(n, to, a, p; t)dt 



M n?=i< 

o(T) 



if M :=n«™« 
otherwise. 



This is the "justification" jCFK + 05l p. 35] for the diagonal term selection, which 
will only keep the terms in (|3"9")l not in o(T) above. 

The expression obtained so far, at the stage (f3"9")h was of the form 

-Pi(B) 



(42) 



A,B s.t. 

AUB=V 
|^|=|B|=fc 



E ( Ott 



2- 



f(A;B)dt, 
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where the crucial feature of / is that it is symmetric in A and B, in the sense that 
/ (ai, • • • ,a fc ;/3i, • • • , (3 k ) = f (a CT (i), • • • , a CT ( fe ); /V(i)> ' ' ' , Ar'(k)) for an Y foe 7 ) G 
x 5fe, if .4 = {a} and 23 = {/?}. In a way that will be made much more explicit 
in Section [51 the diagonal selection step transforms the combinatorial structure of 
(|4"!?]) and thus also (|39[) into an expression of the form 

T fir pl(B) M;S) 



.aub=x> 

|^| = |B|=fe 

where / has the same symmetries as /. Before describing in more details in Section[6] 
the transition from / to /, we will now present the advantages of an expression of 
the form ([43]) . 

4. Laplace Expansions and the Bump-Gamburd Lemma 

Imagine we could obtain the following equality for the numerator in Expres- 
sion ([15]): 

(44) — ) f(A; B) = £ b KiXSK (A) Sx (B), 



i.e. express this numerator as an infinite sum of products of two Schur functions 
s K {A)s\(B), elements of one of the most natural bases for expressions displaying 
this Sk x Sk symmetry. Obviously, the 6 Kj a would still depend on t. We would then 
have, in (|4"3")) . to sum (|4"4")l over balanced disjoint set partitions AUB — T>. However 
we do know that, once this summing is done, we have to recover an expression 
with the full symmetry of (j35|) . i.e. a symmetry associated to permutations of the 
variables of T>, the full group 6>2fc- 

Fortunately, a clever use of the Laplace expansion, due to Bump and Gamburd 
allows to do that very easily. We present a generalized version here. 

Lemma 5 ((see also |BG06j )). Let k and A be partitions, vectorized to length K 
and L respectively. Assume there exists a partition /i that can be vectorized to length 
K + L such thai! 

(45) n + pk+l = sort((K + p K ) U (A + p L )) 

and that \T>\ = K + L. Then, 

s K (A)s x (B) 



(46) w(k,A).* m (Z>)= 



A(A;B) 

AUB=V 
\A\=K 
\B\=L 

where w(/t, A) = ±1 is the sign of the sorting permutation in (|45[) (which is unique 
as p is weakly decreasing). 

If two entries on the RHS of (|45|) are equal, no p can be found as p would have 
to have an increase somewhere. In that case, Equation (141)1 still remains true if its 
LHS is replaced by 0. 



®The symbol U refers here to the concatenation of vectors. As such, (k + px) U (A + p^) ^ 
(A + pi) U («+ Pk)- Once sorted, they are the same, but the sorting permutation will sometimes 
have a different sign. 
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Figure 1. The partitions re = (9, 5, 5, 3, 1) and A = (5, 4, 4, 4, 4, 
2, 2, 1, 1) are (5, 10)-complementary, as the Young diagrams of re 
(represented here in dark blue) and A* = (9, 7, 5, 5, 1) (represented 
here in light red, before its rotation of 180°) can be assembled to 
form a rectangle with 5 rows and 10 columns. 

In other words, provided we can get to an expression such as (|4"4"f. this lemma, 
when used right-to-left, lets us very easily sum over set partitions and recover an 
expression for (|39[) where the full S2k symmetry is obvious. 

of Lemma® It is easiest to start with a partition p, and expand s^(T>) according 
to Equation (jT4j) . The exponents that appear in the matrix on the numerator are 
given by p+px+L- We can use the Laplace expansion on this numerator, separating 
both the K+L rows (= variables) and the K+L columns (= exponents) into subsets 
of size K and L. This gives rise to the partition A U B — T> for the variables and 
the partition of Equation (|45p for the exponents. We thus have the numerators of 
the expansions according to Equation (|T4l) of the Schur functions on the RHS of 
Equation (|46|) , up to a sign. 

The Vandcrmonde in the denominator for s^(T>) can also be split, this time in 
three subproducts: the factors that involve either two variables ending in A, two 
variables ending in B, and one in each. These subproducts end up respectively as 
the denominator in s K (A), S\(B), or, for the "cross terms" in the A(A;B). □ 

This lemma is highly convenient, as we eventually need to take a limit T> — > 
(2k) ■ 0. Since 5^(T>) is homogeneous of degree \p\ in the variables of T>. only one 
term actually matters once we take the limit V — > (2k) ■ of (|35|) . the term where 
fi — 4> is the empty partition. It is thus desirable to know which pairs of re, A will 
produce an empty p in (|45|) . Again, this will admit a very simple combinatorial 
interpretation. 

Lemma 6. Let re and A be two partitions, of partition length less or equal to K 
and L respectively. Then, 

4> + Pk+l = Pk+l = sort((re + p K ) U (A + p L )) 

if and only if A', the transpose partition of A, has vector (L — kk, • • • ,L — Ki). 
Graphically, this means the Young diagrams of k and A* can be assembled into a 
K x L-rectangle. In this case, we also have uj(k,X) = (— 1)' A = (— l) k 

We say that k and A are (K, L)— complementary if they satisfy this last condition, 
which is illustrated by Figure Q] 

of Lemma® This proof follows from a classical combinatorial property |Mac951 
p. 3]: the partitions re and A are (K, L)— complementary if and only if together the 
sequences (m + K — i : 1 < i < K) and (A, + L — i : 1 < i < L) are a permutation 
of the K + L integers between and K + L — 1 . 
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For, the statement about the sign, we first observe that it is true when k = {L K ^ 
is the full rectangle and the partition A is empty. Indeed, in that case no sorting 
is needed. Moreover, if it is true for a pair (k, A), consider a pair of partitions 
A + ), where kT is obtained by removing a box from k and A + by adding the 
corresponding box to A, in such a way that k~ , A + are also (K, L)— complementary. 
In terms of the sequences described above, it is easy to check that this transfer of 
a box merely swaps two of the entries between the sequences. Therefore, it alters 
the sign by 1. By induction on the size of A, we are done. □ 

We will only use complementarity in the case of K = L, i.e. a square. In this 
case, we write Sk(X) for the (k, k)— complementary to A. 



5. Alternate Ending to the CFKRS Recipe 

It will be in our interest to use a slightly different expansion than in Equa- 
tion ([44]) : we will instead expand f(A;B) itself in a symmetric function basis. In 
fact, in Section[6l we will obtain the following proposition, which can be understood 
as a Taylor expansion for /, symmetric in two disjoint sets of variables. 

Proposition 7. Starting from Equation (|39[) . the diagonal term selection in the 
CFKRS recipe leads to Equation (|43[) . with f{A;B) admitting the expansion 
(47) 

f(A;B) = J2 d ^(-l) W s K (A)sx(B) = J2d K xs K (A)s x (B) = cxp f £ W^p^A^iB) 

K,X K,X Kn-u 

where the sums run over partitions k, A or and 
(48) W^:=J2V^ M+H) , 

r>l 

with the VT„ in Q[fc] (defined during the proof, in Section and the F,^ Taylor 

o 

coefficients o/£(s) at r (defined in Section \2.4-ty l- 

Moreover, the d K x are given by polynomials with rational coefficients in the W^ v , 
with total weight in the /i (resp. v) equal to the weight of n (resp. X). 

The proof of this proposition is on page ITT1 

Remark. The constant term of f(A;B), i.e. d^ = exp(W^,p) is well-known, 
and equals the arithmetic factor present in the leading order conjecture of Keating- 
Snaith. We present this derivation in Section 

This last proposition is the missing piece to prove Theorem [3l 

of Theorem^ After its three steps, the CFKRS recipe tells us to evaluate 

(«) f E ( ± y"""%T^^ 

AuB=V 
\A\ = \B\=k 
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We can expand /, thanks to Proposition [7] which takes care of performing the 
diagonal selection, and gives a very explicit expansion for /: 

M k (T)= lim ( T V ^>°^( to g(^))"» 1 ( g ))"^.A^(- 1 )' A| «"(^( g ) di 
aub=v 

\A\ = \B\ = k 

We are now at the crux of the argument in this paper: via iterated applications 
of Pieri's rule (see Section I2.2.4[) , it will be possible and actually easy to express 
S\(B)(pi(B)) u as a linear combination in the basis of Schur functions in the alphabet 
B, i.e. {Sv(B)}. However, thanks to Lemma[5]we know that the outcome of summing 
any pair s K (A)s v (B) over alphabet partitions T> = A U B is actually pretty simple: 
either or a signed s„(D). Actually, thanks to Lemma [BJ we know even more: we 
know which pairs s K (A)s v (B) will produce a ±s M (X>), with n = <f>, which is the only 
/i to really matter once we take the limit T> — > (2k) ■ 0: we need v = This 
can only occur if \v\ = k 2 — \k\, which in turn implies that u = fc 2 — \k\— |A|. 

We put all those ideas together, and get: 
(50) 



Remark that (— l) fc ' k 'o;(k, 5fe(«;)) = 1 by LcmmalHl 



6. Explicit Expansion as Symmetric Functions 



□ 



We still need to explicitly compute the transition brought in our expressions by 
selecting the diagonal terms. We thus need to prove Proposition [7] 

of Proposition [?[ We need to go back to the summands in Expression (|39|) , and ap- 
ply the diagonal term selection to obtain f(A; B). We will first express log(/(^4; B)) 
in the basis of doubly-symmetric functions p ll (A)p u (B) and prove the results about 
the W^v Statements about the d K \ will only come at the very end. 

We first need some new notation. Let Q be the sc10 composed of the inverse of 
primes: Q := {1/p : p prime}. It will be implicit that we have the variable Q run 
through Q. This prevents the overloading of the letter "p" , that would soon occur 
otherwise. 

We can now consider the innermost terms of Expression (1391) : 

iie ^ ■ n e ^ = n (s m^o^-em^^ 

a n>l " m>l Q \u>0 v>0 

Q \u>0 

This was the diagonal selection step: for each Q, the powers of Q have to match, up 
to the slight shifts introduced by A and B. We continue from that last expression, 



"it would be interesting to use Q as an alphabet, but it has proved difficult to do that with 
any advantage. 
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and use first a transition rule described in Section 

nfzM^M^) =n 



K u>0 



K.,A 

\|k| = |-M 



( 



(51) 



JJexp 

Q 



E / K aP k (^)Pa(Q-)Q' 



K.,A 

\I«I = |A|>1 



(52) 



exp 



exp 



E /«a£*(q 



){M}ffi<M-4)®0A(B) 



\|«|=|A|: 



E /kaC({|«|}®^(^) 



\|«|=|A|>1 



7 



The f K \ are defined and computed by using the power series for the logarithm on 
the RHS of the first line, which is valid as = 1. They do not depend on k. More 
information about them is available in Section [9.21 

It is clear from here that we should use the power series expansions for the prime 
zeta function given in Section T2.4. 21 The value of |k| will determine where we base 
our expansion. The value of /mm, and <f>^{A) © <j>m(B) = A(B B are needed to 
derive the next lines, together with the expansions in Section f2. 4. 21 



52J) = exp 



exp 



K, A 

\H=|A|>1 
/ 

E f«xJ2 F $PMA)®MB))- E 1o §(t) 

re, A u>0 7S-40B 

\| K | = |A|>1 



1 



A(A;B) 



exp 



E E F ti 



(m.+n) (m + n 



k,A m,n— 

\| K | = |A|>1 



m, n 



J 



exp 



E f^EK 

K,X 

|k| = |A|>1 



fj, U v 



m M ( K )m t ,(A)|^l|' (K )-' M |6|'( A ^ i ^ ) p AI (^)p I ,(S)/ 



A(A;B) 
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if we use the notational shortcut = ( Ai ... ^ ^ ), remember we had assumed 
|.4| = \B\ = k, and define 

(53) V; v := ( H + H ) J2 /«Am,( K )m,(A)fc'(«)+'W-'M-K^). 

V ^ u ^ / k Ahr 

for any pair of partitions fx, v and r a positive integer. For fixed r, fx, v 1 V is a 
polynomial in /c, as we need > > for the m^n) and m v {\) to be 

nonzero. The exponent of k can therefore never be negative. 

We also deduce from this formula that the deg ree of V^, in k is bounded by 
2r-l(fx) -l{v). 

The d K \ arc obtained by expanding the exponential and using the transition 
matrix between Schur and power polynomials (Equation (|16|) ). They are therefore 
algebraic expressions in the , with rational coefficients and with the total weight 
of the partitions involved in the W less or equal to |k| + |A|. □ 



7. Proof of Theorem [4] 

We wish to prove in this section Theorem |U This is of independent interest. 

Olshanski, Regev and Vershik have proved results on the dimensions of skew 
partitions. Their perspective is a little different from ours, but we can still use 
their results to prove the theorem. 

7.1. Frobenius-Schur functions. Frobenius-Schur functions are used to compute 
the dimensions dim(A, /x). We need to introduce the basics from the beautiful theory 
of [ORV03] . 

7.1.1. Super symmetric functions. Define super symmetric power polynomials in two 
alphabets as 

Pr (A;B) ^E^ + nrE^, 

Px(A;B) :=l[p Xl (A;B) 

i 

We keep the same transition matrices for supersymmetric polynomials as for 
standard symmetric polynomials, so this also defines the supersymmetric Schur 
functions S\(A]B) (via the rule given in Equation (|T5j) ). 

7.1.2. Shifted Frobenius coordinates. Given a partition A, we define the two finite 
sequences of nonncgativc integers pi :— A; — i and qi := \\ — i for 1 < i < d (it turns 
out indeed that they both have length d = d(X), given as the number of entries 
on the main diagonal of the Ferrers diagram of A. This d is called the Frobenius 
dimension of A). This pair of finite sequences is the Frobenius coordinates of A 
noted (pi\qi). For instance, the partition (9,5,5,3,1) of Figure Q] has Frobenius 
coordinates (8, 3, 2|4, 2, 1). The partition (9, 7, 5, 5, 1), in light red in Figure [TJ has 
Frobenius coordinates (8, 5, 2, 1|4, 2, 1, 0). The partition (1) would have Frobenius 
coordinates (0|0) and the empty partition's Frobenius coordinates are (|). 

If A has Frobenius coordinates (pi, ■ ■ ■ ,Pd\qi, • • • ,<ld), define x(A) := (pi + 1/2 : 
1 < i < d) and y(X) := {q t + 1/2 : 1 < i < d). The pair (x(X),y(X)) is then the 
shifted Frobenius coordinates. Note that J2 x ex(\) x + J2 y ey(\) V = W- 
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7.1.3. Definition of Frobenius-Schur functions. Let fj., v be partitions, with m := 
\fi\, n := \i/\. The purpose of the Frobenius-Schur functions is precisely to compute 
dim(A, fi). In fact, this problem serves as their definition, by interpolation of the 
following formula: 

dim(/i, v) = Fs^x^yjv)) 
dim(0, v) n(n — 1) • • • (n — m + 1) 

We caution the reader that Fsa is in general not homogeneous. However it docs 
agree with S\ for its terms of top total degree. 

7.1.4. Giambelli formula. Given two integers a, /, define the hook partition of arm 
length a and leg length I as the partition with Frobenius coordinates (a\l). By abuse 
of notation, we also denote the partition itself by (a\l) then. 

The Giambelli formula, proved for Frobenius-Schur functions in |ORV03[ an 
equation between (2.4) and (2.5)] , is a formula that gives an algebraic expressions for 
Frobenius-Schur functions in terms of Frobenius-Schur functions of hook partitions. 
Let A be a partition with Frobenius dimension d = d(X) and Frobenius coordinates 
(Pii ' ' ' ,PdWi, ■ ■ ■ ,1d)- Then, 

(55) Fs x (A;B) = det [F S(j>iM (A;B)] d id=1 . 

We still need to know how to compute Frobenius-Schur functions of hook parti- 
tions. 

7.1.5. Frobenius-Schur functions for hook partitions. We simply quote |ORV03[ 
Proposition 2.4]: 

P Q 

(56) Fs {pM {A;B) = ^ ^ Cpp'C qq >s {p , W) (A;B), 
where 



(57) 




K V-p' V2'2' ' 2 /' P — Pi 
P' > P- 



The notation t r (A) is a notational shortcut for s^^A) (these functions are also 
part of the elementary symmetric polynomials family). Note that c pp = 1. 

7.2. The proof itself. We are now ready to prove Theorem |4j 

of Theorem^ We know from the results on Frobenius-Schur functions that: 
(5g) dim(«,S fc (A)) _ Fs K (x(S k (X)y,y(S k (X))) 



dim^, S k (X)) (k 2 - |A|) ■ ■ ■ (ft* - \ K \ - |A| + 1) ' 
By symmetry considerations, we also have 

(59) dim(0, S fc (A)) = dim(A, S k (<j>)) = dim(A, (fc fc )), 

as the truncation of the Young lattice to partitions fitting inside a k x k rectangle 
exhibits additional symmetry. 

We use Frobenius-Schur functions again, and obtain 

dim(A, (k k )) _ F 5x (x((k k )); y ((k k ))) 



dim((/), (k k )) k 2 ■ (fc 2 - 1) • • • {k 2 - |A| + 1) ' 
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Therefore, we obtain 
(61) 

f q m , dim(0, (k k ))F SK (x(S k (X)); y(S k (\)))Fs x (x(S k (<l>));y(S k (<l>))) 
dimK5 fe (A)) = fc2. (fea _ 1) ... (fc 2_ N _| A | + 1) ■ 

We will show in the proof of Theorem|4l around Equation (|64|) . that dim(</>, (k k )) = 
g k . Therefore, we are done if we can show the general statement that 

F Sli (x(S k (v));y(S k (v))) 

is a polynomial in k of degree at most 2|/x|. By the Giambelli formula, it is suf- 
ficient to prove this for \i a hook partition. By Equation (|56|). the problem is 
further reduced to supersymmetric Schur polynomials. Finally, by the transition 
rule (fT5j) (and multiplicativity of the power polynomials), we only need to prove 
polynomiality in k for 

p r {x{S k {v)y,y{S k {v))), 

together with a bound on the degree of the polynomial in k of 2r (under the 
hypothesis r > 1). 

It is easy to see that there exist two fixed, finite sets I x and I y such that 

x(S k (v))=I x Ax((k k )) 
y(S k (v))=I y Ay((k k )) 

for all k (with A symmetric difference of sets). Loosely said, up to a few initial 
fixes, the Frobenius coordinates for S k {v) and the full k x k square are the same 
for each k. This implies 

p r {x{S k {v));y{S k {v))) = p r (x((k k ));y((k k ))) - p r (4; /„), 

which is valid for all k. It becomes clear, since the last term in the RHS is a 
constant, that we only need to concern ourselves with the first term in the RHS. 
This term is much simpler, and can actually be explicitly computed: 

p r (x((k k )y iy ((k k ))) = p r (k - \, ■ ■ ■ , i; k - \, ■ ■ ■ , 

= [2£to r is odd; 

1 0, r is even. 

Since Yli=o (*+i) r ^ s a polynomial in k (by Faulhaber's formula) of degree r+1 < 2r 
(since we have r > 1), we are done. □ 

Remark. This proof unfortunately breaks some of the symmetry of the state- 
ment of the theorem, giving a different role to k and A. This is due to the current 
state of the theory on Frobcnius-Schur functions. Careful observation shows that 
we have implicitly proved along the way that 

F 5K (x(S k (X));y(S k (X)))F 5x (x((k k ));y((k k ))) - 
F 5x (x(S k (n)y 7 y(S k (n)))F SK (x((k k ));y((k k ))). 

In fact, one should expect a factorization for the final formula for dim(K, S k (X)), 
that would give a symmetric role to k and A. 
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8. Corollaries of the Main Theorem 

8.1. Leading term. We had actually left out one part of the proof of Theorem 02 
the proof that gk = dim((/), (fc fc )). Indeed, we want our computations to be com- 
pletely independent of previous work of Keating and Snaith for instance. In partic- 
ular, we want to compute the leading coefficient of Pk using only Theorem [3] We 
feel that interesting ideas come from this proof. 

Proposition 8 ((see [CFK + 08] )1. The leading term co(k) of Pk is given by 

afcSfc 



(62) c (fc) 



2M 



with gk also given by 



fc 2 ! 



n-,7= (*+j+ 1 )' 



Proof. For the leading term of Pk , the main sum of Theorem [3] is reduced to just 
one term, corresponding to K — A = <f>. Exploiting that Sk{4>) is the full fc x fc 
square (fc fe ), we have 



(63) c (fc) = -— 



dim(0, (fc fe )) 
Jk 2 



We will show that dim(^>, (fc fe )) = gk and that = a^. 

The combinatorial/RMT factor gk- The new definition given for gk in 
the statement of the proposition is trivially equivalent to the original definition of 
Equation ©. 

The proof that dim(0, (k k )) = j-^-i k ^' + will amount to the hook length 
formula IFRT54] . which gives the dim(</>, A) for a partition A: 

IAI! 



(64) dim(c/>, A) 



with H(X) = Yloex an( ^ the l en gth of the hook based at the box □ of A. 

For a fc x fc square partition, the product of hook lengths is given by 11^7=0 
j + 1), so we are done. 

Equivalently, we are computing here the number of Young tableaux of shape 
fc x fc. This identity between gk and the number of such tableaux was already 
known (see OEIS A039622 [A0311] ). but obscured by the connection with random 
matrices. It is however more natural, as it shows the integrality property of the gk, 
and does not require an extraneous limit in N, as in the definition of Equation ([3]). 

The arithmetic factor <Xfc. We want to show 

k 2 

(65) IT(^ 1_ ~) 2jFi ( fc ' fc; 1; ~) =: ak = ex P(W<M>) = d 4»t> 



The first equality is a reminder of the definition of ak- The last equality is a 
consequence of Equation (|47[) . We thus only need to prove the remaining equality, 
which amounts to logdfc = W$$. 
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We start with Equation (|48|). and obtain 



= E E f^ kl(K)+lw ^ (0) ( b y E q n - <E3» 

r>l Vre.AI-r / 

( 66 ) =/(D(i) fc2 E^r lo gc(")+EE lo r E /KAp K (fc-i) P A(fc-i) 

n>2 r>2 Q \ K,Ahr 

( 67 ) = - k2 EEv ^ -<?")+ E E K E • l )p& • J) 

Q ri>2 71 g r>2 \ K.Ahr 

- fc2 E v ) _ fc2Q + E I Q r E /^«(* • • i) 

n>2 y r>l \ K.Ahr 

(69) =E^ 2 (-Ev) +log fE ^ k ■ • 1)0" 

n>l J \u>0 

Expression (|66| follows from the definition of the Fr an d the definition of each 



p Ki (fc • 1). Expression (|67| rearranges absolutely convergent sums. Expression |68|) 
follows from Mobius inversion, and adds/removes the term —k 2 Q. We go back 
through the definition of the f K \ to obtain Expression (fM)) . 

We are now done, since f)„(fe • I) = ("k^ 1 )- ^ 

Remark. Conceptually, this proof is quite simple: we want the constant 
coefficient in A and B, which means that we will replace A (rcsp. B) by the alphabet 
|^4| • 0, satisfying the properties p K (|-4| ■ 0) = when |k| > 0, p0(|^4| ■ 0) = 1 and 
Po(|-4| ' 5) = I -A I (these values have already been used to obtain Equation (|53l) ~). 
We simply have to be careful (just as CFKRS), since some of the infinite products 
might be divergent if handled inappropriately. In the proof, each of the sums over Q 
converges absolutely, as none of the summands have a remaining Q° or Q 1 term (in 
Expressions ([55)) and there is cancellation in the Q 1 term and no Q° term). 

8.2. Uniform convergence of coefficients. In [RHIlj . Rubinstein and Hiary 
provide uniform asymptotics for roughly the first k coefficients of Pk ■ This relies on 
intricate and separate bounds for the arithmetic and combinatorial components. 

Our method does not bring in significant information for the arithmetic factor. 

On the other hand, the direct combinatorial interpretation we provide in this 
paper might be useful to obtain bounds on the combinatorial factor. In particular, 
Equation {55}, which reduces up to a polynomial the computation of dim(«, Sk(X)) 
to the computation of dim(0, Sk (A)), is remincscent of Lemma 3.3 of |RH11| . which 
reduces up to a polynomial the computation of their pk(ai, ■ ■ ■ ,otu', &k+i, • ■ ■ , &2k) 
to the computation of their pk(X;0). They then translate the computation of the 
Pfe(A; 0) into computation of the polynomials N®(\), via a recurrence formula. On 
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our side, we do not need such a recurrence formula, as the theory of Frobenius-Schur 
functions provide more direct ways to compute the dim((/), S*fc(A)). Indeed, the way 
we presented in Equation (|55|) only requires computing adxd determinant, where 
d = d(X) is the Frobenius dimension (a quantity that grows at most like -\/[A[)- 

This is significant for the results of [RHllj . as the main theorem there only 
concerns the roughly first k coefficients: for a fixed a < 1, the asymptotics in k for 
c r (k) are uniform in the range < r < k a . These asymptotics depend on three 
bounds in [RHllj : Equations (50), (51) and (52). Equation (50) is the only one 
with this type of restriction on the range of r. The other two involve the arithmetic 
factor. It our point of view, it is easily understood why such a restriction might 
appear: when |k| + |A| < fc, we know that dim(K, iSfe(A)) will clearly be nonzero, 
for instance, as k C Sk(X) then. When |k| + |A| ~ k a for a > 1, the situation 
becomes more tricky as the partitions k and A* fitted inside a k x k square might 
start overlapping. 



9. (Pre)computation of the Constants 

The structure of Theorem [3] is such that it relies on several families of constants: 
the Fr , f K \, V^v, W^ v , d^, and dim(K, A). We explain how to compute the 
first five families in this section, as we have already explained in Section [7] how to 
compute dim(«:, A). 

9.1. Precomputation of the Fr . These were defined in (|3~Tj) . pagelTTl 

Their computation only relies on Equation (|32|) , page 1121 which is well-known 
and old. The computation itself involves no new technique, but does not seem to 
appear explicitly in the literature. 

By taking successive formal derivations of Equation (|31[) using the quotient rule, 
we can compute F^ N ' as a series over the squarefree integers q of a (fixed) rational 
function in the Q n '(qr), with n < N (if N equals 0, we actually have to take 
log(£(qr)) instead). If r = 1, when q = 1, the multiple 1 • 1 has to be "corrected" 
according to Equation (|33|) to remove the singularity. 

The main point here is that ^ n > (qr) will be easily computable by Euler-Maclaurin 
summation, and that the terms in all the series involved are decreasing exponen- 
tially fast to 0, so their truncations provide very good approximations. Additional 
tricks can also used to improve convergence, as explained in |Cohj . 



9.2. Precomputation of the f K \. It is clear from their definition in Equation (f5Tj) 
that the f K \ present the symmetry f K \ = f\ K , with the additional restriction that 

/ reA = 0if |/c|^|A|. 

We present some actual values in Table [TJ These were computed using sage 



S + ll] to perform computations in the lazy power series ring over the tensor product 
of the symmetric polynomials algebra A with itself. Sage is helpful as it includes 
all the information necessary to transition between bases of A and therefore also 
A(A) <g) A(B). 

9.3. Precomputation of the VL,. The VL, were defined in Equation (|53|) . There 
is nothing to say actually, except that their computation depends on the previously 
computed f K \, and simply consists of a finite sum. In addition, we have the relation 
V£ v — V£„, which follows from the same symmetry for the f^ v . 
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K, A 


□ 


m B 


= EP ™ 


□ 


i 






m 

B 




i i 

4 4 

1 1 

4 4 




BP 
1 1 1 1 






ill 

9 6 18 

1 1 
6 6 

111 
18 6 9 



Table 1. The values of f K \, with empty entries of value 0, for 
small partition sizes. 



rrm 



B 



mm [J 



D 



1 

16 


1 

12 


1 

32 


1 

16 


1 

96 


1 

12 





1 

24 


1 

12 


1 

24 


1 

32 


1 

24 


1 

64 


1 

32 


5 
192 


1 

16 


1 

12 


1 

32 


1 

16 


11 

96 


1 

96 


1 

24 


5 
192 


11 
06 


11 

192 



Table 2. The values of f K \, for partitions of size 4. 



9.4. Precomputation of the W^ u . We simply truncate the Equation (|48[) . pagefTBl 

(r) 

which defined the W^ v . As this involves the F„ for larger and larger n, this con- 
verges exponentially fast. We still have the symmetry = W vtl . 

9.5. Alternative summation for W M „. Careful consideration of Equation (|4"5)l 
shows that W^ v is obtained as a sum with coefficients of a fixed derivative of the 
prime zeta function, evaluated at successive integers. Sums of this form can be 
efficiently reorganized, leading to new series with much better convergence. This is 
explained in jCohj . close in spirit to parts of [CFK + 08] , 

9.6. Precomputation of the d^ v . These were defined in Equation (|4"T|) . page [TBI 
A similar strategy as for the f K \ works, except that we cannot do the computations 
exactly anymore: the coefficients depend on the W^, which can only be approxi- 
mated. We thus need to change the base field of A (or more accurately A® A) from 
Q to a real field with large precision. Again, the symmetry = d vil is preserved. 
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TP 

n 



EJ 



in 



I 


I 


1 


1 


1 


1 


1 


25 


20 


30 


30 


40 


60 


600 


1 





1 


1 





1 


1 


20 


24 


24 


24 


120 


1 

30 


1 

24 








1 

48 


1 

24 


1 

80 


1 


1 








1 


1 


7 


30 


24 


16 


24 


240 


1 





1 


1 





1 


3 


40 


48 


16 


48 


80 


1 


1 


1 


1 


1 


1 


19 


60 


24 


24 


24 


48 


12 


240 


1 


1 


1 


7 


3 


19 


19 


600 


120 


80 


240 


SO 


240 


600 



Table 3. The values of f K x, for partitions of size 5. 





(6) 


(5,1) 


(4,2) 


(4, I 2 ) 


(3 2 ) 


(3,2,1) 


(3, I 3 ) 


(23) 


(2 2 ,1 2 ) 


(2, I 4 ) 


(I 6 ) 


(6) 


I 

36 


i 

30 


1 

48 


I 

48 


I 

108 


I 

36 


I 

108 


I 

288 


I 

96 


1 

288 


I 

4320 


(5,1) 


1 

30 





1 

40 


1 

40 


1 

90 





1 

45 


1 

240 


1 

80 


1 

80 


1 

720 


(4,2) 


1 

48 


1 

40 








1 

144 





1 

72 


1 

192 


1 

64 


1 

64 


7 

2880 


(4, 1 2 ) 


1 

48 


1 

40 








1 

144 


1 

24 


1 

36 


1 

192 


1 

64 


5 
192 


17 

2880 


(3 2 ) 


1 


1 


1 


1 


1 


1 


1 


1 


1 


7 


19 


108 


90 


144 


144 


324 


108 


324 


864 


96 


864 


12960 


(3,2,1) 


1 

36 








1 

24 


1 

108 


1 

36 


1 

108 


1 

144 





7 
144 


1 

54 


(3, I 3 ) 


1 

108 


1 

45 


1 

72 


1 

36 


1 

324 


1 

108 


1 

324 


1 

108 


1 

16 


1 

54 


131 

6480 


(23) 


1 

288 


1 

240 


1 

192 


1 

192 


1 

864 


1 

144 


1 

108 


1 

576 


1 

192 


1 

144 


17 

4320 


(2M 2 ) 


1 

96 


1 

80 


1 

64 


1 

64 


1 

96 





1 

16 


1 

192 


1 

64 





19 

480 


(2, I 4 ) 


1 

288 


1 

80 


1 

64 


5 
192 


7 

864 


7 

144 


1 

54 


1 

144 





49 
576 


473 
8640 


(I 6 ) 


1 


1 


7 


17 


19 


1 


131 


17 


19 


473 


473 


4320 


720 


2880 


2880 


12960 


54 


6480 


4320 


480 


8640 


25920 



Table 4. The values of f K \, for partitions of size 6. 
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Appendix A. Main Formulas and Definitions 

We now list some of the main formulas and definitions, together with an indi- 
cation of the location where it appeared. It is hoped that this will help the reader 
understand the structure of the final result. 

The symbol "=" indicates a result. The symbol ":=" indicates an implicit defi- 

nition, and the "direction" of this definition. The symbol ":=" indicates an implicit 
definition, provided Conjecture [2] is true. 

Prime zeta function: 

(See 422), p. LH C(s) = 5I^logC("s) 5Rs>l 

n>l 



(See (JSU, p. Q3]) l(r + s) =:J2 F r n) s n reN\{0,l} 

n>0 

(See 431, p.LLU) C(l + s) =: - log(s) + J2 F i" )s ™ 

n>0 

Transition between bases/changes of alphabets/partition invariants: 
(See 48), p.EJ z x := J] (^1%^^) 

(See (HSJ, p.© s A = J2 ~X X (^ 

ph|A| ^ 

(See 4H]), p.© Pai =^ x A (a^a 

Ah|p 

(See 422, p. [TJ3 PmttM) = E ^W^M-WI^"^ 



Final result: 

(See @, p.E) 



(i 



2fc T 



(See 43, p.[3j flfe(x) =: c (fc)x fc + ci(fc)x ~ x + • • • + c fc2 (fc) 

(See@,p.H Cjv(fc) = - ^ d KA dim(A,S fc (K)) 

|k| + |A|=JV 

Dimension: 

(See 47]), p. H 

dim(A, W ) = fc2 - W - I) , gfc g fc) N _ |A| + ir B(*)€Q(*) 
(Expected-to-be-)transcendental coefficients: 

(See 63, p.rjni) X) ^"(^^ :=exp fe^M*)^)) 
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(See (gll), p. ED WV := £ V^^+IH 



r>l 



(See p. E]) VJ,:= 



^ U ^ ' K,AHr 



(See dSJ), p.HHD 



re, A 
|«| = |A 



exp 



\|re|=|A|>l 
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